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Abstract. We consider stochastic model based on the linear stochastic differential equation with the 
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INTRODUCTION 

The presence of l/f noise is ubiquitous in a variety of different systems. Mostly 
l/f noise is Gaussian, but sometimes the signals exhibiting l/f fluctuations are non- 
Gaussian. The non-Gaussianity is often taken as a signature of fluctuator's interaction 
[1]. Nevertheless, statistically independent and noninteracting fluctuators may exhibit 
non-Gaussian noise, as well [2], especially when the fluctuations are strong [3, 4]. 

Recently we proposed stochastic models of 1/f noise based on the nonlinear stochas- 
tic differential equations [5, 6]. The models generate signals with the power-law distri- 
butions of the signal intensity and the power-law spectral densities. 

Moreover, l/f noise is often represented as a sum of independent Lorentzian spectra 
resulting from uncorrelated components of the signal with a wide-range distribution of 
the relaxation times [7]. It should be noted that summation of the spectra is allowed 
only if the processes with different relaxation times are isolated from each another 
[8]. Distribution densities of the signal components described by the linear stochastic 
differential equation are Gaussian and the distribution density of the signal resulting 
from the similarly distributed components is usually Gaussian as well [3]. However, 
the signal consisting of the sequence of components with very different variances may 
be non-Gaussian. In this paper we will analyze the non-Gaussianity of the signals 
exhibiting l/f noise and generated by the linear stochastic differential equations with 
the fluctuating relaxation rate. 

THE MODEL 

Consider the random process x described by the stochastic differential equation 



dx = —y{t)xdt + odW 



(1) 



with the time dependent relaxation rate y(t). Here W is the Wiener process, dW = t, (t)dt, 
with £(t) being the 5-correlated white noise, (^(t)^(t')) = 8{t — t'), and a is the 
intensity (standard deviation) of the white noise. In this paper the stochastic differential 
equations we understand in Ito interpretation. 

When the relaxation rate changes very slowly, we have the signal as a sequence of 
signals with different relaxation rates. 

We can eliminate the parameter o by the appropriate change of the time scale, 
t —¥ o 2 t, while the change of the relaxation rate y(t) may be described by another 
stochastic differential equation, resulting, e.g., in the power-law distribution of 7. 

Therefore, we have the system of two equations, 

dy=a r 'fdW r , (2) 

dx= -yxdt + dW. (3) 

Here oy determines the speed of the change of the relaxation rate 7 driven by the 
white noise and the factor 7" imposes the power-law distribution, P r (y) ~ 7^ (with 
r\ = — 2/i), of the relaxation rate. The diffusion-like motion of 7 should be restricted in 
some interval, e.g., (0,1). Then 

p r (r) = (i+ri)r 71 . (4) 

We can restrict the analysis of the positive x, only, with the reflection of x at x = 0. 
For the definite 7 the distribution of x is Gaussian and the power spectrum is Lorentzian, 

P l (y\x) = 2J^e-^\ (5) 

Si(7l/Wr 2 + « 2 r 1 - (6) 

For very slow evolution of the relaxation rate 7, the resulting characteristics of the system 
(2) and (3) yields from the average of expressions (5) and (6) over distribution (4), 

2(l + r] 



>(x) = J P x {y\x)P r {y)dy 



(7) 



y/nx 3 + 2 v 

S(f) = jPi(y\f)P r (y)dy~i/f-\ (8) 

Here T(a,x) is the incomplete gamma function. 

Therefore, the simple linear stochastic equation (3) with the additive noise, linear 
relaxation and the power-law distribution near zero of slowly changing relaxation rate 
results asymptotically in the power-law distribution of the signal, 

P(x)~l/x 3+2r l, x>l, (9) 

and power-law distribution (8) of the low frequency spectrum. 

For the uniform distribution of the relaxation rate 7, i.e., for /1 = r\ = 0, we have from 
Eqs. (7) and (8) 



1 „ 1 



P(x) = Arerfx ^exp(— x 2 ), (10) 

x L 

S(/)~arctan(l/27T/)//. (11) 





NUMERICAL ANALYSIS 

We have performed numerical analysis of the model (2) and (3), as well. In figure 1 
the simulation results for /i = 0, i.e., for the case of pure l/f noise are presented. We 
see good agreement with the analytical expressions in large intervals of frequency and 
distribution of the signal. Figure 2 shows the dependences of distribution and the slope of 
the spectral density on the parameter ji. In all cases the distribution density of the signal 
exhibits the "fat tail" distributions in contrast to the short-range Gaussian distribution 
for the fixed relaxation rate y. 

We analysed the multifractality of the signals, as well. For this purpose we calculated 
a generalized qih. order height-height correlation function (GHCF) F q {t) defined as [9] 

F q (t) = (\I(t' + t)-I(t')\^, (12) 

where the angular brackets denote the time average. The GHCF F q {t) characterizes 
the correlation properties of the signal /(/), and for a multiaffine signal a power-law 




behavior, 

F q (t)~t H «, (13) 

is expected. Here H q is the generalized qth order Hurst exponent. If H q is independent on 
q, a single scaling exponent H q is involved, and the signal I(t) is said to be monofractal 
[9, 10]. If H q depends on q, the signal is considered to be multifractal. 

Calculation results shown in Fig. 3 indicate that the signal is monofractal with the 
Hurst exponent H w and the slope of the spectrum /3 = 2H + 1, as for random walk. 

CONCLUSION 

The linear stochastic differential equation with the slowly fluctuating relaxation may 
generate monofractal signals with the non-Gaussian 1/fP noise. 
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